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04 ' Abstract. We present high angular resolution images of the H2 1-0 S(l) line emission obtained with the Son of 

^ I ISAAC (SOFI) at the New Technology Telescope (NTT) of the Horsehead nebula. These observations are analysed 

in combination with Ha line emission, aromatic dust, CO and dust continuum emissions. The Horsehead nebula 
illuminated by the 09.5V star a Ori {x ~ 60) presents a typical photodissociation region (PDR) viewed nearly 
edge-on and offers an ideal opportunity to study the gas density structure of a PDR. 
^ ■ The H2 fluorescent emission observations reveal extremely sharp and bright filaments associated with the illumi- 

1^-^ ' nated edge of the nebula which spatially coincides with the aromatic dust emission. Analysis of the H2 fluorescent 

, emission, sensitive to both the far-UV radiation field and the gas density, in conjunction with the aromatic dust 

' and Ha line emission, brings new constraints on the illumination conditions and the gas density in the outer PDR 

region. Furthermore, combination of this data with millimeter observations of CO and dust continuum emission 
allows us to trace the penetration of the far-UV radiation field into the cloud and probe the gas density structure 
Q ' throughout the PDR. From comparison with PDR model calculations, we find that i) the gas density follows a 

^ , steep gradient at the cloud edge, with a scale length of 0.02 pc (or 10") and uh ^ 10^ and 10^ cm~"^ in the H2 

C/3 ■ emitting and inner cold molecular layers respectively, and ii) this density gradient model is essentially a constant 

' pressure model, with P ~4 10® K cm~^. The constraints derived here on the gas density proflle are important for 

^ , the study of physical and chemical processes in PDRs and provide new insight into the evolution of interstellar 

T '~j ' clouds. Also, this work shows the strong influence of the density structure on the PDR spatial stratiflcation and 

- - ' Key words. ISM: clouds - ISM: dust, extinction - atomic processes - molecular processes - radiative transfer - 

horsehead 



1. Introduction 

Photo-dissociation regions at the edge of molecular clouds 
illuminated by bright stars are unique targets for the study 
of the interaction of stellar ultraviolet radiation with dense 
gas. Such interactions ionize and dissociate molecules and 
heat the gas. This is important for the evolution of the 
cloud as well as for the star formation associated with the 
cloud. 

The thermal and chemical structure of the PDRs de- 
pends mainly on two parameters, namely, the intensity of 
the incident far-ultraviolet (FUV, 6< hi' <13.6 eV) ra- 
diation field, characterised by a scaling factor x 

^ The scaling factor x represents here the intensity of the 
incident FUV radiation field in units of the average interstellar 



the hydrogen nucleus gas density, nn- PDR models show 
that due to the attenuation of the incident FUV radiation 
field, the gas temperature and chemical composition are 
spatially stratified between the illuminated edge and the 
inner layer of the cloud. Studying observationally this spa- 
tial structure allows one to derive the local physical con- 
ditions (i.e., X ^-iid i^h) and infer the dominant processes 
in PDRs. Unfortunately for most known PDRs, it is dif- 
ficult to establish the geometry and to precisely measure 



radiation field defined by Draine iDrainel ^T^). For a dis- 
cussion of the definitions of mean interstellar ra diation field 
used in the literature on PDRs see Appendix B of I Allen et alJ 
112004') . 
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the FUV penetration depths and the gas density profile, 
which is a key input parameter of the PDR models. 

In this paper, we study the structure of the Horsehead 
nebula photodissociation region located at the western 
edge of the molecular cloud L1630 illuminated by the 
09.5V star a Ori. This nearby PDR (d-400 pc) viewed 
edge-on offers an ideal opportunity of studying the struc- 
ture of a PDR. It represents one of the sharpest mid- 
infrared (IR) filament (width: 10" or 0.02 pc) detected in 
our Galaxy by the ISO (I nfrared Space Observatory) cam- 
era l|Abergel et allEooij) . In order to resolve properly the 
PDR outer structure, we have obtained imaging observa- 
tions of the II2 1-0 S(l) line emission at 2.12 /im using 
SOFI on the European Southern Observatory NTT tele- 
scope. This line emission is very sensitive to both the FUV 
radiation field and the gas density and the angular reso- 
lution of this observations (1") is ~5-10 times better than 
previous infrare d or molecular and a tomic lines conducted 
in this region JZhou et all Il993t iKramer et alJ. Il996[ 
Abergel et all l2003t IPound et all l2003t iTevssier et al. , 
20041) . Thus, unprecedented constraints on the physical 
conditions at small scale of the illuminated material can 
be derived. Also, in order to trace the PDR spatial strati- 
fication, we compare this near-IR data with the aromatic 
dust and Ha line emission as well as millimeter data of CO 
and dust continuum emission recently taken at the IRAM 
(Instit ut de Radio-Astronomie Millimetrique) 3 0m tele- 
scope l|Abergel et al.U2003l Irbvssier et al.U2004l) and the 
IRAM Plateau de Bure Interferometer (PdBI, IPetv et al.L 
l2004j) . A detailed comparison of the observations with 
PDR model calculations allows us to infer the gas den- 
sity structure throughout the Horsehead nebula PDR. 

The paper is organised as follows. In Sect. 13 we present 
the observations. The PDR models we used are described 
in Sect. 13. ll The results of the calculations which explore 
the effect of varying gas density and of projection effects 
are presented and compared to the observations in Sect. 
I3.2ll3.3ll3.4l and f3.5l In Sect. 01 we discuss our results. Our 
paper conclusions are summarized in Sect.[S| 

2. Observations 

2.1. Horsehead nebula 

In the Ha emission line, the Horsehead nebula (B33) 
emerges from the edge of L1630 as a dark cloud in the 
near side of the HII region IC 434 (see Fig.^). a Ori, a 
09.5V system with Tp ff ~33,000 K. is iii this region the 
main heating source ((Abergel et allEo03h . 

Emissions of the dust and gas associated with 
the Horsehead nebula have been mapped from the 
mid-IR to in il limete r wavelengths llAbergel et all 12005 ; 
Pound et all l2003t iTevssier et all l2004t IPetvetal. , 
2004|) . In the ISO camera (ISOCAM) image taken in 
the LW2 filter (5-8.5 /im), which is dominated by the 
emission of aromatic dust (hereafter referred to as PAHs 
or Polycyclic Aromatic Hydrocarbons), the edge of the 
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Fig. 1. (a) Map of the Horsehead nebula in Ha line 
emission obtained at the 0.9 m KPNO telescope (by Bo 
Reipurth and J. Bally), (b) Map in the L W2 filter (5-8.5 
um) of ISOCAM (6" of angular resolution. lAberael et all . 
\2OOj). The brightness profiles along the cuts marked here 
are shown in Fig. The cuts are rotated by 14° to 
be in direction of the exciting a Ori star situated at 
a2ooo =05h38mni5. 00s, 82000 =-0S" 36'00". 

Horsehead nebula is delineated by a bright and narrow fil- 
ament (see Fig.QJj). The PAH emission, which is found to 
scale with the FUV radiation field (for x ranging between 
1 and 10"^. iBoulanger et all 119983}) . is a direct tracer 
of t he local intensity of the incident radiation. Based on 
this. lAbergel et all l)2003) propose that this filament is due 
to dense material illuminated edge-on by a Ori, so that 
matter beyond the border should not be affected by the 
incident radiation field. Moreover, this structure must be 
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seen nearly edge-on. A large inclination (^30° or more) 
relative to the line of sight would produce a broader fila- 
ment. From the projected distance between the star and 
the interface, dproj ~3.5 pc, and the spectral type of cr Ori, 
the FUV intensity of the incident radiation field is x =60 
(in units of the Draine field, 2.7 10"'^ erg s~^ cm"^). 

The gas density derived from the penetration de pth of 
the F UV radiation is a few times 10'* cm""^ ( Abcrgcl et al.L 
This value is consistent with the analysis (using 
Large Velocity Gradients calculations) of CO emission 
lines obtained at the IRAM 30-m telescope which sug- 
gests that beyond the edge the g as density is of the or der 
of (or superior) to 1-4 lO" cm'^ l|Abergel et alJ.l2003|) . 

However, the angular resolution of ISOCAM (6") is 
still insufficient to resolve properly the density structure at 
the edge of the cloud. Furthermore, the infrared emission 
of the aromatic dust depends on their abundance which 
can yarv throughout the in terface ijBoulanger et al.Lll988l 
Il99(i iBernard et all Il99^ . In this context, we have ob- 
tained high resolution images of the H2 fluorescent emis- 
sion which represents a powerful tool for studying the 
structure of the illuminated gas at the PDR edge. 

2.2. SOFI/NTT observations and data processing 

Using the SOFI/NTT infrared imagcr-spectrometer, we 
have obtained images in the 1-0 S(l) emission line of H2 
at 2.12 fim and in the adjacent continuum filters over a 
section of the western edge of L1630. The image scale is 
0.29" per pixel, giving an observed field of '^5'x5'. We ob- 
served three fields : one centered on the Horsehead nebula 
(^2000 = 5h40mn56s, (52ooo = -2°27'04.9") and two oth- 
ers taken respectively at 110"S 210"E and 125"N 130"E 
of the Horsehead nebula field. The Horsehead nebula field 
was observed in the two band filters H2 1-0 S(l) (Ac=2.124 
lim; AA=0.028 iim) and Br7 (Ac=2.167 ^im; AA=0.028 
/Ltm) over two nights in March 2000, while the two other 
fields were observed over two nights in January 2002 in 
the H2 1-0 S(l), NB (Ac=2.09 fim; AA=0.02 /xm) and Br7 
(only for one field) filters. The seeing was measured to be 
-1". 

Observation and data reduction were performed using 
standard procedures for near-IR imaging. A set of frames 
slightly offset in telescope position were obtained. They 
were then realigned and co-added to produce the final 
image. A calibration was obtained with observations of 
standard stars from the NICMOS hst (N9115, N9116 and 
N9118). Calibration uncertainties are about 10%. 

In Fig. [51 we show a composite map of the three fields 
observed in the 1-0 S(l) H2 line emission with the con- 
tinuum subtracted. To substract from the continuum, we 
have used the NB filter when available, and the Br7 filter 
otherwise. The NB filter is in fact narrower and less likely 
to be contamined by other gas lines. In Fig. O we show 
brightness profiles measured in the H2 1-0 S(l) filter as 
well as in the adjacent continuum filters towards 4 cuts 
(shown in Figs. H and EJ: the cut 1 goes through the main 



bright filament seen at the edge of the Horsehead nebula; 
the cut 2 through the two narrow filaments situated at 
40" N; the cut 3 towards the other bright filament seen in- 
side the L1630 cloud; and the cut 4 through a weak larger 
filament situated inside the L1630 cloud. Concerning the 
accuracy of our measured H2 1-0 S(l) line intensity, we 
note that when using the Br7 filter, the 1-0 S(l) H2 line 
emission subtracted from the continuum can decrease by 
<20%. Hence we assign an accuracy of about 20% to the 
measurements of the 1-0 S(l) line intensity. 

2.3. H2 fluorescent, Ha and PAH emissions 

The SOFI data reveals extremely sharp filaments 
(width:^5" or 0.01 pc) spatially coincident with the aro- 
matic dust emission (see Figs. ^ El and OJ- Since the H2 
fluorescent emission is very sensitive to both the FUV ra- 
diation field and the gas density, this confirms that the 
mid-IR filament seen by ISOCAM at the Horsehead neb- 
ula edge is due to illuminated dense material seen nearly 
edge-on. Towards more tenous material, the H2 emission 
would be weaker and more extended. Also, a large inclina- 
tion relative to the line of sight would produce broader H2 
emission filaments. Furthermore, as expected for a PDR 
seen edge-on, we find that at the edge of the nebula, the 
H2 fluorescent emission peak is displaced further into the 
cloud than the Ha line emission which traces the ioniza- 
tion front (see Fig. Opanel 1). Nevertheless, we note that 
the Ha line emission decreases going towards the cloud 
by only a factor of 5 (see Fig. panel 1). This indicates 
that a Ori is probably located slightly closer to us than 
the Horsehead (if a Ori were beyond the Horsehead, the 
Ha line emission towards the cloud would diminish by a 
much large factor). Moreover, it must be emphasized here 
that since the Ha line emission is a strong function of the 
extinction along the line of sight, it is impossible to make 
qualitative predictions. In cases such as the cuts 3 and 4 
taken towards cloud, the Ha line emission profile probably 
has nothing to do with the appearance of the PDR (see 
below) . 

In the SOFI data, we also discover sub-structures unre- 
solved by ISOCAM. At the edge of the Horsehead nebula, 
we can in particular clearly distinguish two narrow fila- 
ments well separated at several places along the interface 
(see Figs. El and O except around (52ooo=-2°28'00" where 
they are likely superposed. This filamentary structure is 
also seen in the visible (see Fig. Eb) presumably due to 
the scattered light. We believe that these infrared fila- 
ments represent the PDR edge at different positions along 
the line of sight. 

Inside the L1630 cloud, the bright filament seen in cut 
3 may also represent illuminated dense material seen edge- 
on. In fact, the H2 emission is shown to delineate a narrow 
filament spatially coincident with the aromatic dust emis- 
sion. Along this filament, we note that the Ha intensity is 
roughly constant and docs not present any front as at the 
edge of the Horsehead nebula. It is possible that this fila- 
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Fig. 2. (a) SOFI map in the 1-0 S(l) H2 line emission with the continuum subtracted. The seeing was about 1". The 
brightness profiles along the cuts marked here are shown in Fig.|21 (b) ESO-VLT composite colour image with contour 
of the 1-0 S(l) H2 line emission (level is 0.5 10~^ erg s~^ cm^^ sr~^). The seeing was about 0.75". 
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Fig. 3. Brightness profiles measured in the H2 1-0 S(l) filter (solid thick lines, 10^^ erg s^^ cm^^ sr^^) and in the 
adjacent continuum filters - Br7 (dotted lines shown in cuts 1, 2 and 3) and NB (dashed lines, shown in cuts 3 and 
4) - along the cuts shown in Fig. [3 We also show the ISOCAM emission in the LW2 filter (solid thin lines, MJy/sr 
xO.l) and the Ha line emission (dotted-dashed lines, scaled to the H2 emission of cut 1). The offset origins of the cuts 
correspond to the H2 emission peaks. 



ment structure is seen beyond the Horshead nebula. The 
structure seen in cut 4 may be on the other hand due to 
more tenons material since the H2 emission in that case is 
weaker and more extended. Also, we note that the H2 1-0 
S(l) line to LW2 emission ratio towards cut 3 is roughly 
similar to that measured in the cuts 1 and 2 while it is 



lower for the cut 4 by a factor of 2-3 (see Fig. O. This 
may reflect differing local gas densities. The dependence 
of the H2 to aromatic dust emission ratio with gas density 
is discussed in Sect. 13 
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In the next section, in order to study the density struc- 
ture throughout the Horsehead nebula PDR, we combine 
this data with recent miUimeter wave observations. 



(a) Kkm^s 



a* - 




^4l"'l^ Its' ^ O" 40^55' so' 



Fig. 4. (a) Integrated intensity map of the ^"^00(3-2) line 
emission made with the CSO (grey scale) with contours 
of the C^^ 0(2-1) line emission (levels are 4-5, 6, 7.5, 9 K 
km/s) made with the IRAM-30m telescope, (b) Map of the 
1.2 mm dust continuum emission. The brightness profile 
along the cut 1 is shown in Fig. 



2.4. Combination with millimeter data of CO and dust 
continuum emission 

Recently, IPetv et all (|20^) have obtained high spa- 
tial resolution maps of CO and isotopomers emission 
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Fig. 5. Emission profiles throughout the edge of the 
Horsehead nebula along the cut 1 (shown in Figs. QJ 
and^. Upper panel: Ha line emission (dotted lines) scaled 
to the H2 1-0 S(l) line emission (solid thick lines, 10~^ 
erg s~^ cm~'^ sr~^); ISOCAM emission in the LW 2 filter 
(solid thin line, MJy/sr x0.14j sub stracted from the back- 
ground emission taken to be the average of the minimum 
value of this extending cut (6.2 MJy/sr); dust continuum 
emission at 1.2 mm (dashed line, mjy/beamx 0.05). Lower 
panel: Emission of the ^^CO (1-0) (solid lines, K km/s) 
and C^^O (2-1) (dotted lines, K km/sx5) lines obtained 
from with the PdBI. The origin offset is taken where the 
Ha intensity is decreased by 50%. 

lines, namely 12CO(1-0), i2C0(2-l) and 01^0(2-1) at the 
Plateau de Bure Interferometer (PdBI). These data were 
obtained with typical resolutions of 5" and 2.5" at 2.6 mm 
and 1.3 mm respectively for ^^CO and of 6" at 1.4 mm for 
C^^O. We used these data to trace the spatial distribution 
of the molecular gas beyond the edge. 

Secondly, in order to measure the intensity ratios of 
several CO emission lines, we used maps of the {J = 1 — 0) 
and (J = 2 - 1) transitions of ^^CO, ^'^ CO and C^^O 
obtain e d at the IRAM 30-m telescope by lAbergel et alJ 
l|2flfl3|) : iTevssier et all l|2nn4j) . These maps have an angu- 
lar resolutions of 10.5" and 22" at 1.3 mm and 2.6 mm 
respectively. We also use the map of the ( J = 3 — 2) tran- 
sition of CO and i^CO obtained at the CSO (M. Gerin, 
private communication). In Fig. 2Ja), we show the inte- 
grated intensity map of the CO (3-2) line emission with 
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contours of the C^^ 0(2-1) line emission. The comparison 
between these two maps show particularly the effect of 
optical depth in the molecular cloud. The emission of the 
rare and probably opticaUy thin isotopomer C^*0 probe 
the inner shielded regions, while the '^^CO optically thick 
emission likely comes from outer layers. The ^^CO opti- 
cally thick emission could trace the front face-on side of 
the Horsehead nebula. Then, in the following, in order to 
probe the structure in the inner shielded regions, we will 
focus more on the analysis on the optically thin C^^O lines 
emission. 

Finally, in order to trace the quantity of matter to- 
wards the PDR, we used a map of the 1.2 mm dust con- 
tinuum en iission taken w i th MA MBO at the IRAM 30-m 
telescope l|Tevssier et al.U2004) . This map shown in Fig. 
^b) has a resolution of ~11" . We note that the C^^O emis- 
sion coincides generally well with the emission of the dust 
at 1.2 mm, at least in the PDR region where no depletion 
is expected. That confirms that the C^^O probes the inner 
denser regions and that the emission of this isotopomer is 
optically thin. 

In Fig. we show the brightness profiles of different 
tracers across the PDR (cut 1). A layered appearance is 
clear; moving away from the excitation source, the Ha, 
H2, ^^CO and C^*0 lines emissions appear in succession. 
This spatial stratification is typical for an edge-on PDR. In 
front of the PDR ionization front (as traced by the Ha line 
emission), the H2 and CO molecules are photodissociated; 
at optical depth of about Ay < 1, the dissociating photons 
have been sufficiently attenuated to cause neutral H° to 
be transformed into H2; similarly, at Ay ~ 2 — 4, the 
carbon balance shifts from C+ to neutral C° and CO. 
Hence, the separation between the H2 and CO emission 
peaks measures directly the FUV penetration size scale. 
From Fig. [SI we estimate that the observed scale size of the 
region is about '^5-10" per magnitude of visual extinction. 
Further, adopting a hydrogen column density, Nh, per 
magnitude of visual extinction, Nh/Ay = 1.85 10^^ cm~^ 
and a distance of 400 pc, we derive an average gas density 
of ~5 10^ cm~^ for a homogeneous region viewed edge- 
on. This v alue is comparable to the previous gas density 
derived bv lAbergel et all l)2003|) . 

In the following, in order to analyse in details the ob- 
served spatial stratification and constrain the gas density 
structure, we compare the observations with PDR model 
calculations. 



3. Comparison with PDR models 

3.1. PDR models 

We use an u pdated version of the st ationary PDR model 
des cribed in iLe Bourlot et alJ l(l99,'^ (for a recent paper 
see iLe Petit et al.l 120 02*) . In this model, a PDR is repre- 
sented by a semi-infinite plane-parallel slab with an in- 
cident radiation field. The input parameters are (i) the 
incident FUV field X: a-nd (ii) the hydrogen nucleus gas 
density nn. In our models, in addition to the ever existing 



isotropic interstellar radiation field (x=l), we consider a 
blackbody incident spectrum with T^ff — 33,000 K and 
an FUV radiation strength of x=60 (in units of Draine 
field, see Sect. l2.l|l . This in principle is an upper limit and 
hence the effect of changing the value of x is discussed in 
Sect. 01 For the gas density, we will consider constant and 
gradient density models as well as isobaric models. With 
these inputs, the model solves the chemical and thermal 
balance starting from the slab edge at each Ay-step in the 
cloud. 

The dominant process governing the chemical balance 
is the photodissociation of H2 , the rate of which is deter- 
mined by both H2 self-shielding and attenuation by dust. 
F or the dust extinction proper ties, we use the analytical fit 
of lFitzpatrick fc Massal l)l988f) to Galactic average extinc- 
tion curve. The radiative transfer in the absorption lines of 
H2 and CO is treated in detail and the in dividual line pro- 
files ar e described with the prescription of lFederman et all 
We adopt an H2 formation rate Rf (cm^ s ^) con- 
stant throughout the PDR and equal to the standard value 
inferred from observations of H2 UV absorption lines in in- 
terste llar diffuse clouds, i.e. Rj—i 10"^^ cm^ s~^ ((Jural 
This assumption is discussed in Sect.^ 

The chemical network involves about 130 chemical 
species and 1100 reactions. For the elemental abundances, 
we take the following standard \&\\ies, : He/H=0.1, C/H = 
1.35 10-"^, 0/H = 3.19 10-^ N/H=7.5 10-^ S/H=1.85 
10^6^ Fe/H=1.5 10-8. As ^^CO and C^^O observations 
are reported, we have added to the reaction set, the 
main isotopic molecules involving ^'^C (with ^^C/H=1.6 
10-6) and (with i^o/H = 6 10"^) and introduced 
the corresponding fractionation reactions l|Graedel et all 
Il982i^ . Also, the photodissociation of the CO isotopomers 
is treated in details. 

For the physical conditions prevailing in the Horsehead 
nebula (x — 60, uh > 10^ cm^"^, see Sect. 12. it . the ther- 
mal balance is mainly determined by the photoelectric ef- 
fect on small dust grains. The photoelectric heati ng rat e 
is derived from the formalism of .Bakes fc Tielensi 

The gas line emission profiles are calculated assuming 
a plane-parallel edge-on geometry. The gas line intensity 
(in erg s^^ cm^^ sr^^) integrated along a given line of 
sight i can be written as: 

^•'W = A ^ '■PDR 

47r 



(1) 



with J^{i) the line emissivity (in erg s^^ cm^'^) extracted 
from the plane-parallel model at the Av{i)AayeTS, (3 the 
escape probability along the line of sight and I pdr (in cm) 
the length of the P DR along the line of sight. F or we 
use the formalism of lTielens fc HollenbachI jl985f) in their 
Appendix B with a turbulent Doppler width 5vd= 2 km 
s^^. This latter assumption is only important for the ^^CO 
(1-0), (2-1) and (3-2) lines which do not greatly affect our 
conclusions. Optically thin lines (such as H2 and C^^O) 
are unaffected by the assumption used to define (3. 
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We compute the aromatic dust emission similarly by 



IpAH{i) 



JpAH{i) 



An 



X I 



PDR 



(2) 



with JpAH the power emitted by P AHs computed with 
the P AH absorption cross-section of IVerstraete fc Legeil 
l|l992l) . We assume the optical properties and the abun- 
dance of aromatic dust to be constant across the PDR. 
The carbon locked up in PAHs (a <10 A) has an 
abundance of [C/H]pah — 5 10^^ inferred from the 
12 / xm emission per hydrogen in typical Galactic cir- 
rus llBoulanger fc Peraull Il988>l and PDR ()Habart et all 
l200lh and from comparison between observations of 
dust galactic emission and extinction with detailed 



both the values of tih- Fig. El shows also that, as for the 
chemistry, the thermal structure is governed by the pen- 
etrating FUV photons. Until Ay ~ 1, the decrease of the 
gas temperature reflects the dependence of the photoelec- 
tric heating with the FUV radiation field intensity. Then, 
in the inner region opaque to the FUV field, the gas is 
heated by gas-grains collisions, exothermal chemical reac- 
tions and cosmic rays, and so that the gas temperature is 
roughly constant. 



and from comparison 

emi ssion 

mod el calculations ijDesert et al.l Il990l: iDwek et all Il997t 
lllfc Draine. 20olh. 

We compute the 1.2 mm dust continuum emission in 
the Rayleigh- Jeans approximation of the Planck function 
as ( Vlotte et al.„ .1998 ) : 



-(3) 



10 mJy/beam 4 lO^i cm-2 20 K 0.005 cm2g 

with NH-i{i) ~ nH2{i) x I PDR the column density of 
H2 along the line of sight i, Tdustij-) the dust tempera- 
ture and ki, the dust opacity per unit mass. The thermal 
dust t emperature is calculated using the analytic expres- 
sion of HoUe nbach et al ] lll99lh and we adopt fc^ constant 
throughout the PDRs and equal t o 0.0035 cm^ (a s 
recently measured in dense clouds. iBianchi et all l200^ . 
This assumption is discussed in Sect. 01 

To calculate the gas and dust emissions, we assume 
that I PDR (the length of the PDR along the line of sight) is 
constant throughout the PDR. Note that this assumption 
is supported by the fact that H2 and PAHs emissions are 
only seen at the edge of the nebula, so that matter beyond 
the edge is not affected by the incident radiation field. A 
round shape would produce H2 or PAH emission further 
into the Horsehead nebula. 



3.2. Model predictions for a uniform gas density 

As a first step, we consider models with a uniform gas 
density with nn —2 lO'* - 2 10^ cm^"^ in the range expected 
in the Horsehead nebula PDR (see Sects. [TTI and . 

In Fig. El we show the model predictions of the ther- 
mal and chemical stratification as a function of depth 
into the PDR. The chemical structure is characterized 
by the H°/H2 transition occuring at Ay ~ 0.1 and 0.01 
for uh =2 10"* and 2 10^ cm"'^ respectively, and the 
C+/C'^/CO transition occuring at Ay ^ 1 — 4 for both val- 
ues of uh. This can be understood as following. For the 
low x/nn ratio adopted here (~0. 01-0. 001), the H°/H2 
transition is mainly driven by the self-shielding in the H2 
lines. Also, for higher gas density H2 self shields more ef- 
ficiently, so that H° transforms into H2 closer to the edge. 
On the contrary, due to the low abundance of CO in the 
outer region, the C+ZC^/CO transition is governed by the 
dust opacity and occurs at the similar optical depth for 




0. 03 



cl(pc) 

0.003 0.0003 



T,„(K) 





Fig. 6. Model predictions for x = 60 and uh = 2 10'^ 
cmT^ (upper panel) and nn = 2 10"'^ cmT^ (lower panel). 
The visual extinction Ay is taken from the PDR outer 
edge. The scale on the upper axis shows the corresponding 
distance (in pc) from the PDR outer edge. Upper panels: 
Gas and dust temperature. Lower panels: Abundances rel- 
ative to total hydrogen. The chemical PDR structure is 
characterized by the if /H2 and /CP/CO transitions. 
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Fig. 7. Normalised emission profiles predicted by the PDR 
model with a constant gas density equal to uh ~ 2 10"* 
cm"^ (dotted lines) and uh — 2 10^ cm^^ (dashed lines), 
and observed along the cut 1 (solid lines with crosses), as a 
function of the distance from the PDR outer edge, d. The 
scale on the upper axis shows the corresponding offset in 
arcsec. The predicted emission profiles have been convolved 
with gaussian beam profiles of FWHM corresponding to the 
angular resolutions of each observations (given in Sect. 
0). For the observations, d = is taken where the Ha 
intensity is decreased by 50%. 

In Fig. 13 we show the normahsed emission profiles 
predicted by the PDR model as a function of the distance 
from the PDR outer edge. As expected, the emission pro- 
files of both the H2 and CO lines depends strongly on the 
gas density. For higher the emission lines are narrower 
and peaks nearer the PDR edge. In Table Q we give at 
the emission peak the intensity of II2 and CO lines as well 
as several intensity ratios. To estimate IpoR (the length 
of the PDR along the line of sight) we use the II2 column 
density derived from the observed 1.2 mm dust contin- 
uum emission peak (at d ~ 0.02 pc or 10") and adopting 
the dust temperature predicted by the PDR model ^ . For 
n/f = 2 lO'^ cm^^, we find that Tdust ~22 K at d 0.02 pc 
(see Fig.El) and using Eq.Hthat Nh^ - 1.5 10^^ cm'^. For 
riff = 2 10^ cm-3, Tdust ~13.5 K at d - 0.02 pc (see Fig. 
0) and Nh^ ^ 3 10^^ cm'^. This implies Ipdr ~ 0.5 pc 
(or 5') for nu = 2 10"* and 0.1 pc (or 1') for 2 10^ cm'^, re- 

^ To derive Ipdr we do not use the H2 fluorescent emission 
data since this emission could depend on projection effects due 
to sub-structures at the edge of the nebula as shown in Fig. |5| 
and discussed in Sect. 13.51 



spectively. The value of Ipdr derived here for uh = 2 10** 
cm~'^ which corresponds to about 2 times the extent of 
the filaments in the plane of the sky (0.24-0.34 pc or 2-3') 
seems unlikely. If we assume that Ipdr is at most equal 
to the extent of the filaments in the plane of the sky, we 
conclude that uh is greater than 3 10^ cm~^. However, 
this estimate is relevant for the gas into the cloud and 
uh can be lower in the H2 emission zone. For the sake 
of simplicity, we adopt for the model with nn = 2 10^ 
cm~'^ the same PDR length as for nn = 2 10^ cm^'^, i.e., 
IpDR = 0.1 pc. 

TableQlshows first that the absolute intensity of the II2 
fluorescent line depends strongly on the gas density. The 
II2 1-0 S(l) line emissivity peak is two order of magnitude 
higher for nn =2 10^ cm""^ than for nn =2 10'' cm~^ . 
This can be understood as follows. At equilibrium between 
formation and dissociation of H2, we have 

Rf UH - Rd{0) X xe-^' x fs{N{H2)) x uh, (4) 

where Rj is the H2 formation rate, i?d(0) ~ 0.1 Ra{0) 
the unshielded photodissociation rate per II2 for x=l and 
with i?a(0) 5 X 10"^" s~^ the photoabsorption rate, 
e"'^'* and fs{N{H2)) the dust extinction and the H2 self- 
shielding factors, respectively. Since the intensity of the 
1-0 S(l) hue ijj^ is proportional to the photoabsorption 
rate, we have ijj^ cx RfUuN^H'^) where N{H'^) is the 
column density of atomic H atoms. Then, for an edge-on 
PDR we have N{H^) = n^o x Ipdr and as n{H^) ^ uh 
at the PDR edge, we thus have oc n^^. Moreover, 
we see that the II2/PAH intensity ratio scales with the 
local gas density. The emission of aromatic dust Iran 
goes roughly as x e""^"* x Nh x [C/H]pah with [C/H]pah 
the abundance of carbon locked up in PAHs. Thus, for an 
edge-on PDR I pah scales with and ijj^/IpAH oc n^. 

The CO line emission increases also with nn', the opti- 
cally thin lines of C^^O are about 5 times higher for nn =2 
10^ cm~^ while the optically thick lines of ^^CO are about 
2 times higher. The C^^O/^^CO (2-1) hne intensity ratio 
which is mainly dependent on the total column density is 3 
times lower for —2 10** cni~^. Finally, we find that the 
CO(2-1)/(1-0) and (3-2)/(l-0) intensity ratios are slightly 
higher for ni/=2 10^ cm~^. This is probably due to the 
higher collision excitation and the slight increase of Tgas 
in the CO emission region for nH='2 10'°^ cm~'^. 

3.3. Comparison with observations 

We now compare these model predictions to the observa- 
tions. In Fig. [TJ we show the observed emission profiles 
along the cut 1 going through the brightest filament. The 
CO emission profiles have been extracted from the PdBI 
maps. In Table Q] for the observed CO emission lines and 
intensity ratios we use the PdBI data when they are avail- 
able and the IRAM-30m data otherwise. For the purpose 
of comparison of the observations with the models, we 
need to determine the position of the PDR edge. The ion- 
isation front is the natural edge of a PDR but in practice 
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Table 1. Intensities and intensity ratios of H2 and CO lines at the emission peak. 





model 1 


model 2 


model 3 


model 4 


model 5 


Observations 


dpeak (pc) " 














H2 


5 10^2 





0.003 


0.012 





0.01 




0.08 


0.005 


0.016 


0.022 


0.01 


0.02 




0.2 


0.025 


0.03 


0.035 


0.025 


0.023 


Intensity (erg s~^ cm~^ sr~^) ' 














H2 1-0 S(l) 


1 10"* 


1 10"^ 


1 10"* 


6 10"^ 


4 10"^ 


2.4±0.5 lO"'^ 


H2 1-0 S(1)/PAH 


0.01 


0.1 


0.01 


0.008 


0.005 


0.004±0.0015 " 


"CO(2-l) 


4 10"^ 


8 lO"'^ 


7 10"^ 


6 10~^ 


2 10"'' 


8±0.05 10"^ 


(2-l)/(l-0) 


8 


12 


12 


12 


12 


13±1 


(3-2)/(l-0) 


20 


30 


30 


30 


30 


20±2 


"CO(2-l) 


7 10"* 


2 10"'^ 


2 10"^ 


2 10"^ 


2.8 10"'^ 


3±0.06 10"^ 


(2-l)/(l-0) 


7 


8 


8 


8 


10 


11±2 


(3-2)/(l-0) 


15 


20 


20 


20 


20 


13±0.7 


C^**0(2-l) 


2.6 10"* 


1.4 10"'^ 


1.4 10"^ 


1.2 10"^ 


1.5 lO"'^ 


1±0.05 10"^ 


(2-l)/(l-0) 


9 


10 


10 


10 


10 


20±2 



Models 1 and 2 with an uniform gas density equal to nu ~ 2 10* and 2 10^ cm"'^, respectively. 
Models 3 and 4 with a gas density gradient with (3—1 and 4, respectively (see Sect. 13.4"^ . 

Model 5 with a gas density gradient with l3 = 1 and taking into account inclination effect with 9 — 6° (see Sect. 13.^ . 
° Distance of the emission peak from the PDR outer edge. 

' Intensity at t he emission peaks. Fro m the brightness in the LW2 filter (5-8.5 /im), which is dominated by the aromatic 
dust emission jBoulanger et all Il998bl) . we have estimated the aromatic dust emission using the following relationship : 
/pah(2 — 15^m) ~ 2 X 1/1^(5 — 8.5fim) based on ISOCAM-CVF spectrum (corrected from the dust continuum emission). 



marking it is hard. We choose d—0 where the Ha intensity 
is decreased by 50%. 

From comparison with observations, we find that : 

- the H2 emission is over-estimated and peaks too close 
to the edge for both values of njj. This discrepancy is much 
more important for — 2 10^ cm~^ which produces 
a higher H2 line intensity by three orders of magnitude 
and a much narrower H2 filament width. Also, this model 
produces the H2/PAH intensity ratio much higher by a 
factor of the order of 20. For riH—'i 10^ cm"'^, the H2 line 
intensity is over-estimated by a much lower factor (about 
4) and the H2/PAH intensity ratio is roughly comparable 
to the observations. 

- the CO emission peaks much too far inside the cloud 
for nH—2 10** cm"'^. Moreover, for this gas density value 
the C^^O lines emission and the C^^O/i^CO (2-1) hne in- 
tensity ratio are underestimated by a factor of 5 and 2, 
respectively. On the contrary, for — 2 10^ cm"^ the 
peak position of the C^^O emission as well as the CO in- 
tensities and intensity ratios are roughly well reproduced. 

In summary, we find that models with a constant gas 
density equal to riH = 2 10^ - 2 10^ cm~^ cannot re- 
produce the global spatial stratification observed as well 
as the intensities of both the H2 and CO lines. One ex- 
planation of the discrepancies is that the Horsehead neb- 
ula PDR presents at the edge a gas density gradient with 
UH ~ lO'' and 10^ cm"'^ in the H2 emitting and inner cold 
molecular layers, respectively. In the next section, we in- 



vestigate what models with a gas density gradient profile 
predict. 

3.4. Gas density gradient 

We consider models with a gas density gradient profile 
defined as : 

n^(d) = I < ^ i^i^ (5) 

[n^H : d>dQ ^ ^ 

with d the distance from the edge of the PDR and n'j^ 
the constant gas density in the d > do region. We take 
10^ cm"'^ since high value of nn is required to ac- 
count for the observed C^^O emission profile, as discussed 
above. For the scale length (do) and the slope (/?) of the 
gradient, we adopt do =0.02 pc (or 10") corresponding 
to the peak of the 1.2 mm dust continuum emission and 
/3=l-4. 

In Fig.lHl we show the model predictions of the thermal 
and chemical stratification as a function of the depth into 
the PDR and in Fig.|^the corresponding emission profiles. 
In Tabled we give the corresponding intensities and inten- 
sity ratios. It clearly appears that the models with a gas 
density gradient reproduce much better the observations. 
In fact, for these two models, we have nn ~ lO'' and 10^ 
cm"'^ in the outer and inner region, respectively, and the 
emission profiles of both the H2 and CO lines arc found to 
be roughly similar to the observed ones. Moreover, these 
models predict intensities of both the H2 and CO lines 
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emission, as well as the intensity ratios, comparable to 
the observations. 

For both models considered here, one can see that the 
scale length of the gas density gradient (dg) corresponds 
to that of the gas temperature change. Tgas decreases in 
fact from T""^^ -200-300 K to T™™ ~10-20 K on a scale 
length of the order of —0.02 pc (see Fig. ISJ. At the edge 
where nn is low the gas is —10 times warmer than in 
the inner region where uh is 10 times higher. Also, we 
find that for these models the thermal pressure is roughly 
constant throughout the PDRs and equal to — 4 10^ K 
cm~^. This suggests that the thermal pressure dominates 
the gas pressure. 

In order to test this interpretation, we run a model 
with a uniform pressure equal to 4 10^ K cm~^ as found 
above. In Fig. ^| we show the model predictions of the 
thermal and chemical stratification as a function of depth 
into the PDR. The gas density goes from — 10^ cm~^ in 
the H2 emitting region to 5 10^ cm"'^ in the inner cold 
molecular layers and the thermal and chemical structures 
are very similar to that of the gas density gradient models 
(see Fig. IS)). For the isobaric model, the scale length of 
the gas density and temperature gradients is also of the 
order of do ~ 0.02 pc. Consequently, we find that the 
corresponding emission profiles and intensities are in good 
agreement with the observations as for the gas density 
gradient models. 

Nevertheless, we note that for both gas density gradi- 
ent and isobaric models considered here there are several 
discrepancies with the observations: 

- the H2 emission profile is narrower than the observa- 
tions and the H2 intensity is overestimated by a factor of 
2-4. Moreover, the H2 and PAH emissions have a tendency 
to peak too close the edge. 

- the C^*0 line emission is shifted too much in the 
inner region by < 0.01 pc (or < 5"). 

However, projection effects along the line of sight may 
be important and could affect the width and the position 
of the emission profiles, as well as, the absolute inten- 
sities. In the following, we investigate the impact of the 
projection effects resulting from PDR inclination and sub- 
structures. 



3.5. Projection effects 

First, we discuss the effects of PDR inclination on the 
emission profiles. For a PDR tilted by an angle relative 
to the line of sight, the emission integrated along a given 
line of sight i can be written as: 



cl(pc) 

0.01 0.003 



3=^ 



An 



X l{j)/sin{9) 



(6) 



where j labels the various layer crossed by the line of sight 
(see Fig. Ill|l . J{j) is the line emissivity extracted from 
the plan-parallel model at the Av{j)-\ayeTS, l{j)/sin(9) 
the length of layer j along the line of sight and jmax is 
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Fig. 8. Same as in Fig. a model with a gas density 

gradient with f3 = I (upper panel) and (3 — 4 (lower panel) . 



defined as Ej=r Kj)lsm{e) = Ipdr/cos{0). For the CO 
line emission, J{j) is multiplied by f3{j) the escape prob- 
ability of the layer j along the line of sight. 

To examine the effect of the inclination, we calculate 
the maximum angle of inclination 0max considering that 
the PDR edge presents a wall of density. Using the ob- 
servations of H2, 9max can be given by sin~^{dH2 /IpDii) 
with ^ 0.01 pc the distance of the H2 emission peak 
from the PDR edge. Thus, taking IpoR = 0.1 pc, we find 

In Fig. El we show the results for the gradient density 
model with (3=1 inclined by 6* = 6°. As expected, the gas 
and dust emission profiles are displaced further near the 
edge (by —0.005 pc or 2.5"). Moreover, the H2 emission 
profile is larger (FWHM x 2) and the H2 intensity is re- 
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Fig. 9. Same as Fig.^^for models with a gas density gra- 
dient with (3 — \ (dotted lines) and (3—4 (dashed lines). 
For more details see text in Sect. \3.4\ 
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Fig. 10. Same as in Fig. \^ for a model with a uniform 
pressure equal to 4 10^ K cm~^. 



duced by a factor of 2 (see Table ^1. The PAH emission, 
which is more extended, is on the contrary practically not 
affected and the H2 /PAH emission becomes similar to the 
observed one. Finally, note that for the other models con- 
sidered here the inclination effects are found to be similar. 



These results show that the inclination effects could ex- 
plain the discrepancies between the observations and the 
gas density gradient and/or isobaric models. Nevertheless, 
taking into account these effects the H2 and PAH emis- 
sions peak still too close to the edge. This discrepancy 
could be in part due to the projection effects resulting 
from the filamentary sub-structures observed at the edge 
of the Horsehead nebula (see Sect. 12. 3() . In fact, for the 
comparison with the model, we have selected a cut across 
the interface where the two filaments appear closest but 
probably they are not precisly aligned along the line of 
sight. Furthermore, the position of the PDR edge is not 
really well determined and it also has a slant. 




Fig. 11. Schematic diagram of a PDR inclined by an angle 
9 relative to the line of sight. Index j lables the Ay -layers 
of the model. 



4. Discussion 

In summary, we conclude that in order to explain the ob- 
served spatial stratification of the different tracers, as well 
as their intensities and intensity ratios, a high gas density 
gradient at the edge of the Horsehead nebula is required. 
The scale length is about do ~0.02 pc (or 10") and the 
value of uh is about 10'' and 10^ cm^'^ in the H2 emit- 
ting and inner cold molecular layers, respectively. Since 
the scale length of the gradient corresponds to that of the 
gas temperature change (Tg^s decrases from ~200-300 to 
10-20 K on ^0.02 pc) and the thermal pressure is roughly 
constant throughout the PDRs, we conclude that the gas 
density is given by the thermal pressure. A model with a 
uniform pressure equal to ^4 10^ K cm^'^ fits in fact also 
well the observations. 
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Fig. 12. Same as Fig. ^ for models with a gas density 
gradient with (3 — 1 (dotted lines) and taking into account 
the inclination effects with an angle 9 — 6° (dashed lines). 

This study is consistent with the earher findings of 
such a gas density gradient on the western edge of the p 
Ophiu chi molecular cloud bv lHabar'3 l)200l|) : lHabart etaP 
l)2003j) . For this interface (x ~250), we find from the anal- 
ysis of the observed profiles of the H2 1-0 S(l), C° ^Pi-^Pq 
and ^■^CO(3-2) emission lines that the gas density profiles 
is as in Eq.[Slwith n^^ ~4 10'' cm'^, p ~2.5 and do -0.05 
pc corresponding to the scale length of the gas tempera- 
ture change. 

Here it must be emphasized that our results are based 
on several important assumptions. 

i) We assume x=60 while x could be lower. However, a 
large angle between the sky plane and the line connecting 
the Horsehead to the star seems unlikely (see Sect. l2.ll and 
I2.3|l and, even diminishing the value of x by a factor of 2 
will not cause significant differences. 

ii) We assume that the H2 formation rate, Rf, is equal 
to the standard rate derived in the diffuse ISM while it 
could be higher. In Ha,bart c t al. (200j), we show that 
for moderately excited PDRs (i.e., x <1000), such as the 
Horsehead nebula, the H2 formation rate is a factor of ~5 
times larger than the standard rate, R^. However, we have 
run models with i?j = 5 x i?° and find that for the condi- 
tions prevailing in the Horsehead nebula the chemical and 
thermal structure is practically not affected. Nevertheless, 
note that the H2 1-0 S(l) line intensity will be increased. 
However, as our results are mainly based on the analysis 
of the spatial stratification, we conclude that the assumed 



value for Rf does not greatly affect our conclusions. 
Hi) We adopt k^, constant and equal to 0.0035 cm^ g^^ 
which has an unc ertainty of roughly a factor 2 
ijBianchi et all EoO^fl . However, a change in the value of 
ki, by a factor of 2 will only modify IpoB. and the line 
intensities by at most a factor 2. 

iv) There are also two fundamental uncertainties in our 
present study. One is that we assum e a homogeneous gas 
density while PDRs could be clumpy* Aberg el et alJ lj2003t) 
conclude that the material behind the illuminated edge 
may be non homoge nous with clump sizes of — 10""^ pc 
or smaller. However, iGorti fc HollenbachI l)2002|) have re- 
cently investigated the effects of FUV radiation on clumps 
and show that in PDRs, the ambient FUV penetrates 
through the surface of a dense clump and heats this sur- 
face to high temperatures, causing mass loss and thereby 
inducing photoevaporation of the clump. Smaller clumps 
are shorter lived compared to larger clumps and the pho- 
toevaporation timescale would be as short as 5 10'^ yr 
for a clump size of 10~^ pc and an average mass of ~ 10~^ 
M0 (or a density of 10^ cm~'^). 

v) Another assumption is that of a static, equilibrium 
PDR. In reality, the propagation of the ionization and 
photodissociation fronts will bring fresh H2 into the zone 
emitting line radiation. The effects of a moving ioniza- 
tion front become very important for the PDR struc- 
ture if the timescale for H2 to flow across the PDR is 
smaller than the timescale for H2 photodissociation, lead- 
ing to the merging of the dissociation front with the 
ionization front. For a column density of H2 from the 
PDR edge of N{H2) ~ lO^^ cm'^ (which is typically 
our case), the criteria for merged the dissociation fron t 
to the ionization front is (Ho llenbach fc TielensL Il99fl|) : 
^IF ^ 3 (x/i^h) km s^^ with vip the velocity of the 
ionization front. Apparently, the condition is not met in 
the Horsehead nebula, since the observed layered struc- 
ture shows that the H2 offset from the ionization front 
(see Fig.O. Qualitatively, the merger tends in fact to oc- 
cur for PDRs heated by hotter stars with a high r atio o f 
extreme-UV to FUV photons l|Bertoldi fc Drain! Il99^ . 
Nevertheless, it is possible that the dissociation front is 
shifted closer to the ionization front beca use o f the advec- 
tion of H2. However, .Storzer fc HollenbachI l)l998() have 
modelled the PDR structure assuming an ionization front 
moving into the PDR and found that non equilibrium ef- 
fects are probably minor in objects similar to the Orion 
Bar, where the ionization front propagation was expected 
to have significant impacts. We conclude that presently, 
it is reasonable to use stationary PDR models in order to 
estimate the physical conditions. 

5. Conclusion 

The aim of this study has been to determine the gas den- 
sity structure of the Horsehead nebula PDR viewed edge- 
on using the observed emission profiles of different tracers. 
We have obtained high angular resolution (~1") imaging 
observations of the H2 1-0 S(l) line emission, using SOFI 
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at the NTT, in order to resolve the PDR structure in 
the outer zone. Analysis of the H2 fluorescent emission, 
very sensitive to both the FUV radiation field and the gas 
density, in conjunction with the aromatic dust and the 
Ha line emission, brings strong constraints on the illumi- 
nation condition and the gas density in the outer PDR 
region. Moreover, in order to trace the penetration of the 
FUV radiation field into the cloud and probe the gas den- 
sity structure throughout the PDR. we have combined the 
infrared emission of H2 and aromatic dust with millime- 
ter observations of CO and dust continuum emission (ob- 
tained at the PdBI and the IRAM-30 m telescope). 

From comparison with PDR model calculations, we 
show that a constant-pressure PDR model with P ~4 10^ 
K cm~^ best fits the global spatial stratification observed, 
as well as, the absolute intensities and intensity ratios of 
the different molecular lines and dust emissions. Constant 
density models fail to reproduce all the available data at 
the same time. A high gas density gradient with a scale 
length of about 10" (or 0.02 pc) and njj rising from ~ 10^ 
cm~^ in the H2 emitting region to ^ 10^ cm~^ in the 
shielded cloud interior is required to explain the observa- 
tions. The gradient scale length corresponds to that of the 
gas temperature change {Tgas decrases from ^200-300 to 
10-20 K on -^0.02 pc); the gas density gradient model cor- 
responds to a constant-pressure model and we conclude 
that the PDR is in pressure equilibrium (at P ~4 10^ K 
cm~^). 

To the best of our knowledge, this study represents 

the first observational evidence of the dominant infiuence 
of the thermal pressure on the gas density structure of 
a PDR. Experimental evidences of the thermal pressure 
domination arc extremely scarce up to now because of the 
small spatial scales involved. The constraints derived here 
on the gas density profile are important for the study of 
physical and chemical processes in PDRs and provide new 
insight into the evolution of interstellar clouds. 
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